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It is shown that an effective theory with meron degrees of freedom produces confinement in SU(2) 
Yang Mills theory. This effective theory is compatible with center symmetry. When the scale is set 
by the string tension, the action density and topological susceptibility are similar to those arising 
in lattice QCD. 
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I. INTRODUCTION 

Although confinement is one of the most striking and 
fundamental phenomena that arises from the deceptively 
simple QCD Lagrangian, its physical mechanism has yet 
to be understood. One important element is sufficient 
disorder to drive a large Wilson loop to produce an area 
law as seen, for example, in the strong coupling expan- 
sion of lattice QCD, which produces an area law already 
in lowest order. However, disorder alone is not the whole 
story, as evidenced by the fact that the strong coupling 
expansion also erroneously yields confinement in U(l) 
gauge theory. A second important feature is center sym- 
metry. By gauge fixing in continuum QCD up to a resid- 
ual center symmetry or in lattice QCD by multiplying 
all the links in the time direction on a single time slice 
by an element of the center, it follows that the effective 
action has the symmetry S e s(P) — S c s(ZP) for an ele- 
ment of the center, A serious candidate for the 
confinement mechanism should have the potential of real- 
izing this symmetry. One appealing analytical approach 
to understanding nonperturbative QCD is expansion of 
the path integral for the partition function around sta- 
tionary classical solutions and evaluating the fluctuations 
around these solutions. Development of highly success- 
ful instanton liquid models and observation of instantons 
and their zero modes in lattice QCD have provided clear 
insight into how chiral symmetry breaking arises in QCD 
UQI. However, singular gauge instantons fail to produce 
confinement and, as will be seen below, also fail to pro- 
duce center symmetry. 

In this work, we reexamine the possibility that merons, 
another set of solutions to the classical field equations, 
can produce confinement. As noted long ago [5| , merons 
are sufficiently disordering that they have the potential 
to produce confinement. In contrast to instantons, whose 
gauge fields fall off as r -3 in singular gauge, merons can 
only be written in regular gauge with fields that decrease 
as r _1 . In the absence of correlations between distant 
merons, this long-range gauge field would give rise to an 
unphysical background field. Since analytic treatment of 
these correlations appears to be intractable, the behavior 
of meron ensembles has not previously been analyzed. 

In SU(A) gauge theories, the deconfinement transition 
is associated with a symmetry property. Center symme- 



try, a residual discrete gauge symmetry, is realized in 
the confined phase and broken in the deconfined phase. 
In SU(2), center symmetry transformations change the 
sign of the Polyakov loop, which we define here along the 
compact z-direction of Euclidean space, 

P(t, x, y) = ^ tr P exp{i <j> dzA 3 (t, x, y, z)}. 

Realization of the center symmetry implies vanishing 
of the expectation value of the Polyakov loop which in 
turn guarantees an exponential decrease in the Polyakov 
loop correlator related to confinement. In the deconfined 
phase, the Polyakov loop expectation value is finite. The 
Polyakov loop thus serves as an order parameter. In 
the construction of ensembles of field configurations with 
which to explore confinement, it is useful to study the 
Polyakov loop of the building blocks. As Fig. 2] shows, 
unlike the winding number, the asymptotic value of the 
Polyakov loop of an instanton is not changed in the tun- 
neling process and is identical for instantons and anti- 
instantons. Thus, tunneling processes do not give rise to 
a coherent superposition of fields with opposite values of 
P so that a center-symmetric ensemble of field configu- 
rations cannot be generated in a natural way by super- 
position of singular gauge instantons. In contrast, the 
tunneling process described by a single meron changes 
the sign of the Polyakov loop as can be seen from the 
figure. Furthermore, the asymptotic values of P for a 
meron-antimeron pair and a meron-meron pair are of op- 
posite sign. Thus, merons do not single out either one of 
the two center-elements of SU(2) and have the potential 
of generating a center symmetric ensemble. 



II. MERON EFFECTIVE THEORY 

In this work, we explore the idea that merons are the 
essential degrees of freedom in SU(2) Yang-Mills gauge 
theory by writing the partition function as a path integral 
of an effective action depending on the positions and color 
orientations of an ensemble of merons, 



Z = I dzjdhj 



T S[A(zi,hi)] 



(1) 



The gauge field for a meron in Lorentz gauge with its 
center at the origin, after appropriate choice of the coor- 
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FIG. 1: Polyakov loop P(t, 0,0) as a function of time. The 
solid line corresponds to one instanton centered at the origin. 
The dashed and dotted-dashed lines show the results for a 
meron-meron pair located on the f-axis at t = — 7, —2 and an 
antimeron-meron pair at t = 2, 7. 

dinate system in color space and after regularization of 
the singularity, is given by 



(2) 



Color and space-time dependence are correlated via the 
't Hooft tensor r\ alll , 0. For vanishing meron size p, a(x) 
is a solution of the Euclidean classical field equations [6| • 
Antimcrons differ in sign if one of the space-time indices 
is 0. The action density of merons and antimerons is 
given by 



s ( x ) — 2 tr -f)"^f - ,.J x r -2 , ; 



[.x 4 + 4xV + 8p 4 ]. 

(3) 

Unlike instantons, the field-strength decays asymptoti- 
cally as 1/x 2 , giving rise to an infrared logarithmic sin- 
gularity in the action. For vanishing meron size, the ac- 
tion is also logarithmically divergent in the ultraviolet, so 
we will use p as an ultraviolet regulator. The topological 
charge density of a meron or antimeron 



s(x) = ±~tr F^F^ = ± ^^ [3 



2p 2 



(4) 



leads to a finite, size-independent topological charge 

"=8^2 J ^ i xS(x)=±± (5) 



The meron ensembles to be considered in this study 
contain field configurations obtained by superposition of 
merons and antimerons of fixed and equal number Nm /2, 



N M 



A n( x ) = ^ h^a^x - z(i))h 1 (i 



(6) 



Such a configuration is specified by the position of the 
centers z(i) and the color orientations 

h{i) = ho(i)+ih{i)-cr , h 2 (i) + h 2 (i) = 1. (7) 



In the ensembles to be discussed, the location of the 
merons is restricted to a hypercube 

-1<*„(»)<1, V = 16. 

We identify the effective action in with the Yang Mills 
action. Our standard choice for meron size and coupling 
constant is 

p = 0.16, g 2 = 32. 

It is essential to note that the infrared divergence of the 
action of a single meron does not prevent construction 
of physical ensembles having an extensive action with a 
large number of merons. For instance, with the following 
choice of the color orientations of a system of 4 merons 
(or 4 antimerons) 



h k {i) 



-i,k ■ 



(8) 



the action density decays as 1/x 6 . 

The meron ensembles have been generated by Monte 
Carlo sampling of the action in the path integral, Eq. 
||TJ. In each step of a Metropolis update, the position 
and color orientation of a given meron are tentatively 
changed, the induced changes in the action density are 
evaluated at a set of mesh points distributed over the 
whole volume, and the configuration is accepted or re- 
jected based on the global change in action. The long 
range nature of the meron fields makes the changes ex- 
tend throughout the whole system. 

We first present the results for the central quantity in the 
discussion of confinement, the Wilson loop 



U = -tr jPexpi J dx^A^x) 



(9) 



The integral is ordered along the closed path C. Our stan- 
dard choice is a rectangular path located in an (x,i,Xj)- 
plane with the center at the origin and with the ratio 
of the sides equal to 2. For a given configuration, we 
evaluate twelve different Wilson loops Wij and obtain 
our final results by taking the ensemble average and the 
average over the twelve orientations. The statistical er- 
rors are calculated from the variance of the twelve ori- 
entations. Fig. [2] shows the universal behavior of Wil- 
son loops calculated with meron ensembles with a range 
of meron numbers. After rescaling the area A — > A ,4, 
the values of the Wilson loop lie on a universal scaling 
curve, where deviations from scaling are within the sta- 
tistical uncertainties. For loops that are not too small, 
the universal curve can be parameterized by the sum of 
a perimeter term and an area term 

\n(W) ^uj + tVxV -aXA (10) 
as shown in Fig. [5] The values of the parameters are 



a; = -0.72, t = 0.93, a = 8.17 



(11) 



The concave shape of ln(VF) as a function of the area 
arises from the difference in sign between the area and 
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FIG. 2: Logarithm of Wilson loop as function of area A. The 
data correspond to meron ensembles at g 2 — 32 with Nm = 
1000, 500, 200,100, and 50. The area has been rescaled by 
1.86, 1.41, 1.0, 0.75, and 0.53 respectively. The line shows 
the fit with Eq. iT(l . 
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TABLE I: Meron density, meron size, action density and 
topological susceptibility for meron ensembles containing Nm 
=1000, 500, 200, 100, and 50 merons 

perimeter terms and makes the presence of a positive 
string tension unambiguous. Numerical results also con- 
firm the expected increase of the Wilson loop with in- 
creasing perimeter at fixed area. These results thus 
demonstrate confinement in these meron ensembles. 



III. TOPOLOGICAL SUSCEPTIBILITY AND 
ACTION DENSITY 

In view of this impressive scaling behavior, it is useful 
to scale all subsequent results with a physical scale set 
by the string tension. For convenience in thinking about 
physical magnitudes, although we are considering SU(2), 
we will identify the value of a with the physical value 4.2 
fm -2 of the string tension. Ultimately we will compare 
dimensionless ratios involving powers of a with SU(2) 
lattice results, so that this arbitrary physical scale in no 
way affects our quantitative results. Thus, our unit of 
length (u. 1.) is converted to physical units by 



1 u.l. = 1.4VAfm. 



(12) 



Table 0] summarizes our principal results for the action 
density (s), and the topological susceptibility \ calcu- 
lated in meron ensembles for different values of the meron 



density tim = Nm/V and scaled to the string tension. A 
striking feature is the fact that the action density and 
topological susceptibility depend essentially on a single 
scale, which we have chosen as the string tension. In- 
deed, despite the changes in the number of merons by up 
to a factor of 20, and of meron size and meron density 
by factors of 2 and 3 respectively, the action density and 
topological susceptibility vary by less than ±5%. The re- 
sults turn out to be similarly insensitive to large changes 
in the coupling constant. Hence, we will write the action 
density and topological susceptibility as 



(a) = 13.1 a 2 ~s{np\g 2 ) , 

X V4 = 0.44^/2^1/4(^2) 



(13) 



where, within five percent s « 1, x 1 / 4 ~ 1 . In the regime 
of parameters investigated, breaking the scale invariance 
of the meron fields by introduction of the meron size p 
appears to have an effect similar to that of introducing a 
finite lattice size in the quantum theory. In both cases, 
the underlying scale invariance manifests itself in the ap- 
pearance of logarithmic singularities for vanishing coor- 
dinate space regulators. 

To understand the interplay between meron size and 
meron density we observe that the action density in the 
center of a meron [cf. Eq.Q] is much larger than the 
average action density in the meron ensembles 



s(0) _ 0.87 
Is) 



pV 2 



= 13.5 



(14) 



p=0.25 fm 



This suggests separation of the action density into back- 
ground and meron peak contributions, 



sb + s M ■ 



(15) 



The meron peak contribution is obtained by integrating 
s(x) [Eq. over a sphere of radius r, which for small 
meron size {p r) becomes 



SM 



n / d i xs(x) 



37T 2 7J 



5 

12 



In ■ 



(16) 



The matching requirement on r 

s(r) = s B 

yields the following expression for the action density 
3^ 2 



(17) 



SB 



671" /t> 

I - 

4 V3 



5,2 4 



This expression makes explicit the logarithmic singularity 
of (s) with the meron size in the small p limit if sb is 
identified with the physical (i.e. regularized) value of the 
action density. With the action density s as an input 
(cf . Table HJ) , the relative strength of the background 
contribution can be determined numerically 



-^=0.65 



0.75. 
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This result indicates that a significant fraction of the ac- 
tion density is associated with the logarithmically singu- 
lar contribution in Eq. Ijlfijl . 

We can now compare the meron results in Eq. I|13|) with 
known QCD results. The topological susceptibility 

x= (si) 2 / d M^V^(*) f; v f;m)- (is) 

is a robust quantity that can be measured unambigu- 
ously in our effective theory and in lattice QCD. It has 
a direct physical interpretation because of its relation 
to the rf mass by the Veneziano-Witten formula. The 
most extensive lattice measurement of the SU(2) topo- 
logical susceptibility by extracting the continuum limit 
from calculations over a large range of the coupling con- 
stant yields x 1/4 /o- 1/2 ~ 0.483 ± .006. The meron 
result x 1/4 /c 1/2 ~ 0.44 from Eq.(E3 is in excellent agree- 
ment with this lattice result. Note that x can be com- 
puted reliably for the meron ensembles since the topo- 
logical charge of a single meron, Eq. ©, is finite. The 
topological susceptibility is dominated by the short-range 
peaks in the topological charge density associated with 
individual merons and antimerons, and not by the long 
range background field. The contribution due to single 
peaks yields 

x 1/4 = 0.505 ntf. (19) 

The weak variation in x 1 ^ 4 an d the magnitude of this 
estimate are in qualitative agreement with the results of 
Table |H 

Because of the necessity of subtracting divergent terms 
to define the continuum limit, the action density, or 
equivalently the gluon condensate, is more difficult than 
the topological susceptibility to evaluate accurately in 
either QCD or a meron ensemble. Lattice SU (2) calcula- 
tions for the action density range from (s) /a 2 = 4.5 8] to 
25.3 Bj and QCD sum rule results range from (s) /a 2 = 
4.5 [T(| to 10 [Tj. The meron result (s)/cr 2 ~ 13.1 from 
Eq. i|13|) is thus consistent with our present knowledge of 
the action density. It is useful to note that in the lattice 
calculations of Ref. [g, the divergent contributions are 
about a factor of 20-500 larger than the extracted value 



of (s). As noted above, the meron action density also 
contains a divergent contribution that, if subtracted, will 
significantly reduce the value of the condensate given in 
Tabled 



IV. CONCLUSIONS 

In conclusion, we have shown that an effective theory 
with meron degrees of freedom succeeds in describing es- 
sential features of QCD. The long-range gauge fields pro- 
vide the correlations and disorder needed for confinement 
while the short range fluctuations play an essential role 
in other physical observables. They render the meron dy- 
namics more complicated than for ensembles of weakly 
interacting singular gauge instantons and must therefore 
be treated numerically. The central feature and a major 
success of this approach is producing confinement. We 
have demonstrated a confining area law and discussed 
how a meron ensemble can implement center symmetry. 
We also observe scaling behavior, reminiscent of lattice 
Yang Mills theory which might not a priori have been 
expected. In our effective theory there are three param- 
eters, g , p, and um, and one combination of p and xim 
is determined by fitting to the string tension at a given 
g 2 . The action density and topological charge density 
are rather insensitive to the other combination and the 
Wilson loops scale to a universal curve. This scaling be- 
havior is presumably connected with the scale invariance 
of the multi-meron action in the limit p — > 0. In addi- 
tion, we have shown nearly quantitative agreement with 
the value of the topological susceptibility as measured 
in SU(2) lattice QCD and have obtained a value of the 
action density of the correct order of magnitude. 
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